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Abstract 

We consider super symmetric (SUSY) models in which a very hght gravitino 
is the hghtest SUSY particle. Assuming that a neutralino is the next-to-lightest 
SUSY particle, we present a measurement of the neutralino mass at the LHC in 
two photons -|- missing energy events, which is based on the Mt2 method. It is 
a direct measurement of the neutralino mass itself, independent of other SUSY 
particle masses and patterns of cascade decays before the neutralino is produced. 



Among various supersjTiimetric (SUSY) models, those with an ultrahght gravitino of 



gravitino problems [T] . In this letter, we assume a neutralino is the next-to-lightest SUSY 
particle (NLSP), and present a measurement of its mass at the LHC. It is based on 
the so-called Mt2 method [2]. We show that this method can directly determine the 
neutralino mass, independently of other SUSY particle masses, and it does not rely on 
specific patterns of cascade decays before the neutralino is produced. 

In the scenario considered here, essentially all the SUSY events will end up with two 
neutralino NLSPs|l] each of which then dominantly decays into a gravitino and a photonic 
We assume that the decay length of the NLSP neutralino is so short that the decay 
occurs inside the detector and the photons' momenta are measured well. Therefore, the 
main signature at the LHC will be two high transverse momentum photons and a large 
missing transverse momentum. If such a signal will indeed be discovered, one of the most 
natural candidates for the underlying model is a SUSY model with a gravitino LSP and 
a neutralino NLSP. 

Furthermore, from the prompt decay of the neutralino, we can assume that the grav- 
itino is very light, essentially massless for the following discussion. This is because the 
NLSP decay length is proportional to the gravitino mass squared as 



and a heavier gravitino {171^/2 > 0{1) keV) would make the neutralino decay outside the 
detector Ifl This indirect information of the massless LSP plays a crucial role in the NLSP 
mass determination. 

Let us start by briefly explaining the Mt2 method [2j. Suppose that there is a particle 
A which promptly decays by the process A B + X, where B is a visible (standard- 
model) particle and X is a neutral and undetected particle. When two A s are produced 
in a collider, we can measure the two i?s' transverse momenta Px'^, Px'^ and the missing 
^We assume _R-parity conservation. 

^We do not discuss the case in which the neutrahno mainly decays into a Higgs/Z-boson and a 
gravitino. 

^For a moderate gravitino mass corresponding ctnlsp = 0(10) cm — 0(10) m, the neutralino decay 
causes "non-pointing" photons [3]. The present method of the neutralino mass determination may also 
work in this case. 



mass m3/2 < 0(10) eV are very attractive, since they are completely free from notorious 
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transverse momentum p™'^^ = p^'^ + p^'^B The Mt2 variable is then given by 

(Mt2)'= min fmax|(M^'V,(M4'V|l, (2) 

T "r/-^T ~/^T 

where the minimization is taken over all possible momentum splittings, and 

(M^'V = m| + + 2(E™E^'' - ■ p?'') for t = 1,2 (3) 



with E^^^^'^ = yrrix + |p™^^'*P and E^'^ = yrn^ + |p^'*p. This Mt2 variable is designed 
to have the endpoint at itla when we input the correct value of rrix- However, in general, 
the mass mx of the missing particle X is unknown, and therefore one can obtain only a 
relation between mx and rriA^ 

A crucial point in the scenario considered here {{A, B, X} = {neutralino, 7, gravitino}) 
is that we can assume the massless LSP, mx = = 0, as discussed above. Therefore, 
we can directly determine the NLSP mass by the Mt2 method. As we show in the 
Appendix, the Mt2 variable in this case is analytically expressed as 

(M = I "^P^P^^ ci < or C2 < . . 

^ \ for ci > and C2 > ' 

where pi = Ipt^I, P2 = IPt^I; and C2 are given by 

pr^ = cip^^ + c2P^', (5) 

and 2; is a real positive solution of the following equations: 



A{a-by = (^2(a + 6) + 3- l)(J2(a + 6) + 3 + 3) 



\3 



1 

r 3 



f 3 



T 1 

COS 6' + Ci sin^ 9- 

2 z 

cos6' + C2 sin^ 9- 

2r z 



7,1 7,2 

r = — , cos 6'= . (6) 

Pi P1P2 



^We assume that the missing pr^ is dominantly caused by the two Xs and the contribution of other 
sources of missing p^p are neghgible. 

^See also recent developments in the Mt2 method [4]. 

^ For completeness, we also show an analytic expression of Mt2 for the case of massive LSP {mx ^ 0) 



in the Appendix. 
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Figure 1: Mass spectrum of SIGM 
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Note that Mt2 is completely defined by the missing transverse momentum and photon 
momenta, independently of other kinematical variables. We should also emphasize that 
the present method does not rely on a direct pair-production of the NLSPs, i.e., we do 
not assume back-to-back transverse momenta of the NLSPs, p™'^^ -|- p^^ + pij.'^ = 0. In 
the following, we show how this method works at the LHC, by taking explicit examples 
of gauge mediated SUSY breaking (GMSB) models which realize the mass spectrum with 
an ultralight gravitino LSP and a neutralino NLSP. 

We consider two gauge mediation models for a demonstration. In the following, mass 
spectrums are calculated by ISAJET 7.72 [5j and we use programs Herwig 6.5 ^ and 
AcerDET-1.0 [7J to simulate LHC signatures. 

The first example is a strongly interacting gauge mediation (SIGM) model [8J, in 
which the NLSP is a neutralino. We take the same SIGM parameters as the example in 
Sec. 4 of Ref. [8]. The mass spectrum is shown in Fig. [H The masses of the lightest 
neutralino and gravitino are 356 GeV and 10 eV, respectively. 

We take the events cuts as follows: 

• > 4 jets with pt > 50 GeV and pt,i,2 > 100 GeV. 

• > 2 photons with > 20 GeV. 
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Figure 2: A distribution of Mt2 for the SIGM example, (a): parton level signature, (b): 
detector level signature. 



• Meff > 500 GeV, where 

Meff = EPT,+P^^ (7) 
jets 

• pf^^ > 0.2Meff. 

Under these cuts, we see that the standard-model backgrounds are almost negligible. 

In Fig. [2]-(a), a parton level distribution of Mt2 is shown for an integrated luminosity 
of 10 fb~^. Here, we take the sum of gravitino and neutrino transverse momenta as the 
parton level missing p-p. As discussed in Ref. [S] , very little number of leptons are produced 
in the SIGM. Therefore, missing is due to almost only gravitinos and the assumption 
that p™'^^ = plf^^ + plf^'^ is satisfied. There is a clear edge at Mt2 — m^o = 356 GeV. 

In Fig. [2]-(b), we show a distribution of Mt2 after taking account of detector effects. 
In order to extract the point of the edge, we use a simple fitting function; 

f{x) = {ax + b)9{~x + M) + {cx + d)9{x - M), (8) 

where 6{x) is the step function and a, b, c, d and M are fitting parameters. We fit the data 
with f{x) over 300 < Mt2 < 500 GeV and find 

m^o = 357 ± 3 GeV. (9) 
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Figure 3: Mass spectrum of SPSS 

Here, the estimation of the error is done by 'eye' because of lack of information on the 
shape of the Mt2 distribution. The estimation that m^o = 357 ± 3 GeV is very good 
agreement with the true value m^o = 356 GeV. 

Next we show another example. We study the Snowmass benchmark point SPSS p], 
which is a minimal gauge mediation model with a neutralino NLSP. In Fig. [31 SPSS mass 
spectrum is shown. The masses of the lightest neutralino and gravitino are 139 GeV and 
4.S eV, respectively. 

In Fig. Ill-(a), a parton level distribution of Mt2 is shown for an integrated luminosity 
of 10 fb^^. The event cuts are the same as in the previous SIGM case. The blue and 
dashed line represents the case that p™'^^ = I]gravitino Pt the red and solid line p™*^*^ = 
Sgravitino Pt + I^ncutrino Pt' SPSS, there are many neutrino production sources. Hence, 
we cannot see a clear edge as in the SIGM case. However, there is a cliff at Mt2 — ttl^o = 
139 GeV. 

In Fig. Il]-(b), detector level distribution of Mt2 is shown. To get the value of m^o, we 
fit the data with f{x) in Eq. (ED over 110 < Mt2 < ISO GeV. Then we get 



m^o = 139 ± 3 GeV. 

Xi 



(10) 
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Figure 4: A distribution of Mt2 for the SPSS, (a): parton level signature. The blue 
and dashed line represents the case that p^^^^ = Z^gravitino Pt ^^e red and solid line 
Pt'"" = EgravitmoPT + EneutrmoPT- i^)- detector level signature. 



The error estimation is done by 'eye'. This value agrees with the true value [m^o = 
139 GeV). 

In summary, we have presented a determination of the neutralino mass for the SUSY 
models with an ultralight gravitino LSP and a neutralino NLSP, which may work in the 
early stage of the LHC 

Though we have considered GMSB models with a neutralino NLSP, our method is 
applicable to any model in which the signal events will lead to a pair of cascade decays 
that result in 

any cascade decay ^ A ^ B + X, (11) 

where i? is a visible (standard-model) particle and X is a missing particle that is almost 
massless. The mass of A is then determined by the two Bs' momenta and the missing 
transverse momentum. 

For example, let us consider GMSB models with a slepton NLSP. In this case, the 
slepton, lepton and gravitino correspond to A, B and X in Eq. (II ip . respectively. In 
addition to leptons from the sleptons' decays, many other leptons are produced in this 
scenario. However, we may see which of observed leptons is produced through the slepton 
decay by measuring lepton's momentum, or by detecting a kink of its track for a long- 
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lived slepton. In such a case, we can measure the slepton mass with the Mt2 method as 
discussed above. Furthermore, the present method may work in an axino LSP scenario. 
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Appendix 

In this appendix we derive Eq. (j4j). We start from Eqs. ([2]) and (l3l). We assume that B 
is massless. 

(i) mx = case: First, we consider the case that X is a massless particle. The Mt2 
variable is defined by Eqs. ([2]) and ([3]) with rriB = ttlx = 0. If a momentum splitting is 

ji 1 • miss,l X,l 1 miss, 2 X.2 iiti • 11 

the correct one, i.e., = Pt and = Pt' ? then each transverse mass is smaller 

than the mass of A, uia'- 

m\ = (p^'^ = 2(|p^'^||p?1 cosh(Ay^) - p^'' ■ p?'^) > (M«)^ (12) 

for i = 1,2, where Ay* is the rapidity difference of B and X in each decay chain. From 
this, it is clear that 

Mt2 < mA. (13) 

We do not assume the relation 

Pt'' + Pt'' = -pT^ (14) 

which holds in the case of a "back-to-back" pair production of A s. We may assume that 
Pt'^ and Pt'^ are linearly independent and p^^^^ can be expressed as 

= cip?'^ + C2p?'i (15) 



8 



Here, ci and C2 are real coefficients and they are given by 



Cl 



sin^ 9 



BA 



B,2 



cos 6 



P1P2 



where 



C2 



sin'^ 



miss B .'2 miss B.\ 



cos 6 



P1P2 



B,l B,2 

I I -8,2 1 r, Pt ' Pt 

Pi = Pt 5 P2 = Pt > cosfc* = . 

P1P2 

The momentum sphtting p™^'''''^ and p™'^^'^ can also be expressed as 



miss.l / \ B,l , B,2 

Pt = (ci - a;)PT + VPt > 



miss, 2 BA , / \ ^ 

Pt = a;PT + (c2 - 1/)Pt 



-B,2 



where x and ?/ are real variables. We rewrite Eq. ([2]) as 

(Mt2)^ = 2pip2 min \max{ zi{x,y), Z2{x,y)} 

x,y&ii. L 1^ J 

where 

{M^'\x,y)f 



zi{x,y) 



2piP2 



^1 ^ ^ , /ll 2 , 9-2/1 

h ?/ COS & +2/ sm 



Cl — X 



y cos 



{M^'\x,y)y 



Z2{x,y) 



X COS 9 + {c2 — y)r]^ + x"^ sin^ 9 — [x cos 6* + (c2 — 



2piP2 

and r = p2/pi. It is clear that 



> 0, and zi{x,y) = <^=^ |/ = 0& x<ci, 



Z2{x,y) > 0, and Z2{x,y) = -v^ x = Sz y < C2. 



From this, we can infer that 



(Mt2)^ = if Cl > & C2 > 0, 



2 fl) 2 f2) 2 

and for other values of ci and C2, (Mt2) is given by {Mj. {x,y)) = (M^ (x,?/)) at the 
point (x, y) = (xo, yo) where the contours of zi{x, y) and -22(2;, y) in the x-y plane become 
tangent to each other. We denote the corresponding value z = zi{xo,yo) = Z2{xQ,yo) in 
the following. Eqs. ([22]) and ([23D yield 

rsin^6'/ zcos^\^ rz 

^0 = ^ [yo-—^] + o • 2/1 27 

22; V sm / 2 sm 



sin 6 / z cos ^ \ ^ z 

2/0 = — o ^0 ~ + o • 2n + C2 (28) 

2rz \ sm 9 J 2r sm 9 



with the tangential condition 

sin'^6'/ 2; cos 6* \/ 2; cos 6* \ , , 

-^(^o-—^)[yo-—^) = i. 29 

z"^ V sm 6* / V sm 9 / 

We can obtain z by solving these three equations. 

A straightforward calculation yields that these equations reduce to 



4(a - by - {^2{a + 6) + 3 - l)i'sj2{a + 6) + 3 + 3)^ = 0, (30) 

where 

1 /r 1 \ 

a = — I cos 9 + Ci sin^ 6^- ) , (31) 



r 3 



1 



2 



6 = rt( — -cos^ + C2sin^^-V (32) 
V2r zj 

It can be checked that the above equations have a unique real positive solution of z. 
Eqs. (l30l) -( |32l) have been used for the analysis in this work. 

In the special case of a "back-to-back" pair production, in which Eq. ( IT4l) holds, we 
recover the result obtained by taking the massless limit of the formula in Ref. [TO] . 

(Mt2)\ ^ ^ ^ =2(|p?'^||p?'Vp?'^-p?'2) = 2pip2(l + cos^). (33) 

back— to— back 
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(ii) nix 7^ case: Generalization of the above result for the case with massive X, i.e., 
rrix 7^ 0, is straightforward. In this case, the Mt2 variable is defined by Eq. ([2]) with 
TUB = 0- The same argument as above shows that Eq. (fT3|) holds also in this case. 
Calculating in the same way as above, it can be shown that 

(Mt2)' =mj, + 2pip2Z (34) 



with z being the solution of Eq. ( 1301) with 



1 /r ^ .9^1 rsin^6'm^v. 1 \ 

— (--cos^ + cism^^ 3^—), (35) 

r 3 



2 z 2 pI z'^ 



6 = r3 — -cos^ + casm^^ ^ • 36 

V2r 2; 2r 2;^/ 

For the case with massive X, this expression for Mt2 is valid for any values of ci and C2. 

The existence of a unique positive real solution of z can also be checked. 
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